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ABSTRACT . With standing wave ultrasonic techniques, small changes In phase velocity which result from 
changes In some external parameter (e.g., temperature or magnetic field) have traditionally been determined by 
observing shifts in the mechanical resonance frequency of a composite resonator. Some previous investiqators 
have assumed that the fractional change in velocity Av/v is equal to the fractional change in frequency Av/v. vie 
discuss quantitatively the errors involved in such an approach. We show that using the relation Av/v = Av/v 
results in substantial, inaccuracies when the loading effect of the transducer(s) cannot be neglected. Substan- 
tially improved formulas for determining the dispersion are presented and one of these is shown to be much more 
accurate than all jDrevious approximations. The results of simulated and actual experiments over wide ranges of 
dispersion Av/v, , transducer loading parameter 6 = JL t and frequency are analyzed in order to compare 

the errors inherent in the various approximations. t c s s 


I. Introduction 


In a variety of ultrasonic experiments one 
monitors changes In acoustic phase velocity in a speci- 
men that result from variations in some external para- 
meter, e.g., magnetic field, pressure, or temperature. 
Standing wave ultrasonic techniques are well suited to 
measurements of this sort. Fractional changes Av/v 
in phase velocity can be determined from measurements 
of the shifts in frequency of a standing wave mechani- 
cal resonance. To analyze the data, the relation 
Av/v a Av/v is conmonly used, where Av/v is the frac- 
tional change In mechanical resonance frequency. 

Because measurements are usually made on composite 
resonators consisting of a specimen plus one or two 
transducers, this "uncorrected formula" is only approxi- 
mate. The uncorrected formula results in substantial 
errors in the estimation of the ultrasonic dispersion 
Av/v for a variety of cases of experimental interest. 

In this paper, we present substantially more 
accurate formulas for determining the dispersion Av/v 
from standing wave ultrasonic measurements. Using nu- 
merical simulations we compare quantitatively the . 
errors resulting from use of the new approximations 
with those resulting from use of the uncorrected formula. 
In order to demonstrate the significance of the new 
dispersion formulas, experimental data are analyzed 
using the conventional and the improved formulas, and 
the resulting values for the dispersion are compared. 

II. Theory 

A. Reflection Case (One Transducer) 

Me consider an ultrasonic resonator consisting 
of a specimen (properties labeled with subscript s) and 
one transducer (subscript t). The velocity of sound 
Vg is given by 

Vs » v”/n , (1) 

where A, is the length of the sample and v” is the 
frequency of the ntn (sample only) mechanical resonance. 
In the limit that the transducer has no effect on the 
mechanical resonance frequencies of the sample, vg is 
equal to v^. the measured mechanical resonance frequen- 
cy. Thus, 

avs/vj " avj/vj + aij/tj . (2) 

Limiting the discussion to cases where the last term in 
Eq. (2) can be neglected, one obtains the "uncorrected 
formula" for the dispersion 

• (3) 


Equations (2) and (3) ignore the loading 
effects of the transducer. For the case where the 
transducer loading parameter 6 = P+At/Pc^s 
large, a more accurate expression for the mechanical 
resonance frequency of the composite resonator is“ 

v", = - vj) . (4) 

where is the (unloaded) transducer resonance fre- 
quency. Using Eqs. (l)-(4) and restricting to the case 
where v" is not too far from v^, one obtains after 
some manipulation the "(1+6) formula" for the disper- 
sion. 


Av^/v^ = (AvJ/v”)(1+6) . (5) 


Anticipating the results of Section III, Eqs. (3) and 
(5) typically exhibit errors substantially larger than 
those resulting from experimental inaccuracies. 

In order to develop a more accurate formula 
for the dispersion Avc/v^, we begin with the one trans- 
ducer resonance condition® 


tan 


27t£._v_ 

s c 


Pt^ 

Ps''s 


( 6 ) 


Here T = tan{Trv^./v.f-) , and and are the densities 
of the transducer and sample, respectively. For con- 
venience we have suppressed the superscript n which 
specifies the particular mechanical resonance which is 
being monitored. After a change in the external para- 
meter (e.g., magnetic field), the same resonance equa- 
tion relates the new velocity v| to the new composite 
resonance frequency vj. 
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where T* = tan(7ivi/vf ) . (We are assuming negligible 
changes in tj.) ^ 

We seek the value of vj, given the values of 
'^s* Ahe remaining parameters. Subtracting 
Eq. (6) from Eq. (7) and simplifying, one obtains 
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PfVi I j » 

where X = — ~ - — j , The argument of the tan- 
gent function on the left side of Eq. (8) is small ovei 
a wide range of parameters, so we keep only the first 
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term In the series expansion. The resulting equation 
Is quadratic In v|. 


A{>^)^ + B V* + C = 0 , 

(9a) 

where A • + Tfiv^H, 

(9b) 

B • -T* v*5v4.M - NT - V vi, 
St sc 

(9c) 

C • NT*v^ 
and where ^ 

(9d) 

W • sin(4irfcgv^/v^) 


> “ cos^(2it£^ 
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(10b) 

The resulting one transducer formulas for vj and the 
dispersion Av^/Vg are given by 

. -B * / B^-4AC 
s ZA 

(11a) 

AVg/Vs * (vj - Vg)/Vg . 

(11b) 


form with Eqs. (9) and (11), but requires the following 
redefinitions of H and N,. replacing Eq. (10): 

M • sin(2Ti£gV^/v^} 

+ I C05^(7Tt^v^/v^) (15a) 

N =» 05ti) 

The dispersion formula defined by Eqs (9), 
(11), and (15) does not apply to the cotangent set of 
resonances. The corresponding formula for detf rminin 
the dispersion from transmission measurements tn a cc- 
tanqent resonance [Eq. (14)] is given by Eqs. (9) and 
(11) with a second redefinition of M and N; 

M « s1n(2irt^v^/Vg) 

+ |sin^(irtgV^v^) (16a) 

N = sin(2Trt^v^/v^) . (16b) 


(The plus sign in Eq. (11a) is required.) 

In Section, III, we investigate the behavior of 
the uncorrected formula [Eq. (3)], the (1+6) formula 
[Eq. (5)], and the present result [Eq. (11)] using both 
nunerlcal simulations and experimental data. 


B. Transmission Case (Two Transducers ) 

For a composite resonator consisting of a 
sample and two transducers, the appropriate resonance 
equation Is® 


where 


tan 85 - 2 tan 8^ 

2 

* (tjI) tan =■ 0. (12) 

e, -2irV^/v, . 


* ”^''t 


r - (PjV, - Ptv^)/(0,v^ * p^v^). 

This equation can be factored into two simpler reso- 
nance equations which are individually quite similar to 
the resonance equation for the one transducer case 
[Eq. (7)]: 


tan 

cot 
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The sol utl OTIS to Eq. (13) (the "tangent set" of reso- 
nances) describe only alternate menfcers of the full set 
of resonances implicitly defined by-Eq. (12). The solu- 
tions of Eq, (14) (the "cotangent set" of resonances) 
represent the other half of the full set of resonances. 

Using Eq. (13) and analytical techniques simi- 
lar to those outlined for the one transducer case, an 
Improved dispersion formula for the tangent set of 
resonances is obtained. The result is identical in 


In order to analyze data in the two transducer 
case, one must determine whether the mechanical reso- 
nance being monitored is a member of the tangent cr the 
cotangent set. One first calculates the number C, 
given by 

C » v^/AVg , (17) 

where Av^ = Let C = [C], the greatest integer 

less than or equal to»C‘. If C is odd, the first 
resonance Vr above is a tangent resonance. If C is 
even, the first resonance above is a cotanaent 
resonance. This allows an unequivocal determination of 
whether dispersion monitored with any particular mechan- 
ical resonance should be calculated usina Eqs. (9), (11) 
and (15), or using Eqs. (9), (11), and (16). In Section 
III we Investigate the behavior of the various aooroxi- 
matlons to over a wide range of conditions. 

III. Discussi on 


In this section, we examine the behavior of the 
uncorrected formula [Eq. (3)] , the (l+d) formula 
[Eq. (5)] , and the present result [Eqs. (9)-{ll)]. A 
numerical simulation technique is used to fnvestioate 
the errors that result from the use of each of the three 
approximations to AV5/V5. Data from an actual disper- 
sion experiment are also analyzed. Since the dispersion 
formulas for the transmission case [Eqs. (9), (11), and 
(15) or (16)] exhibit a behavior very similar to that 
for the reflection case, only the results of the analy- 
sis of the reflection case [Eqs. (9)-(n)] are presented. 
In order to study analytically the errors re- 
sulting from the various approximations, . we used compu- 
ter iteration to find to 5 parts in 10'’ the solutions 
Vc of Eq. (6) with assumed values for v^, pj, vt, 
p^, and Vf. An assumed value of dispersion Avj/vs = 

(v| - Vjl/Vg defines the shifted phase velocity v^. 
Repeating the iteration - process with this new value 
for the phase velocity yields values for the shifted 
mechanical resonance frequencies v?. The approximate 
formulas [Eqs. (3), (5), and (9)-(ll)] were used to 
obtain estimates of AvjAs for each vc. vj Pair. The 
percent error for each approximation is defined with 
respect to the assumed value of AV5/V5.. . The para- 
meters were chosen so that the errors of the various 



approximations could be studied as functions of the 
magnitude o,f the dispersion AVg/Vg, the size of 6,* and 
the distance in frequency of the mechanical resonance 
Vc from vf. 

In Figures 1 and 2 we present the results of 
our analysis of Eqs. {3), (5), and {9)“(11), with the 
magnitude of the dispersion AVg/Vg as the independent 
variable. Figure 1 treats the case where Vc is close 
to vt- (In this case Vc was chosen to be the first 
mechanical resonance above v^.) Results for two 
values of 5 are presented: 6 small (0.005), a value 

typical of experiments in solids, and 5 large (0.2), a 
value typical of liquid experiments. In any specific 
experiment, 6 is constant, so the three curves for 
each value of 6 are to be viewed as a group. 



Figure 1. Absolute value of the percent error in 
Avs/vs versus AVg/Vg for each of the three approximate 
formulas and for two values of 6 = Pt^t/Ps^s* using 
the first mechanical resonance v^. above v^. 

For 6 large or small, the uncorrected formula 
(Eq. (3)] and the (1+6) formula [Eq. (5)] exhibit 
errors that are essentially constant over a wide range 
of AVj/Vc. This result Is discussed below. For 
6 « O.OOi, use of the uncorrected formula results in 
errors of approximately 0.5%, while use of the (1+6) 
formula results 1n errors of about 0-01%. By contrast, 
the present result [Eqs. (9)-(11)) is far superior, 
yielding errors from 3 to at least 6 orders of magni- 
tude smaller than either of the other approximations. 
Since data for Vq and accurate to parts in 10® or 
10^ are available from experiments, ’ the Increased 
accuracy provided by Eqs. (9-11) is required. (Since 
the final accuracy of our double-precision computer 
calculations is limited to about 10‘®%, we plot a 
dashed line where the error becomes less than this 
value. In this region, we can only establish an upper 
bound for the error.) For the larger value of 6 = 0.2, 
the errors for all three approximations Increase, but 
the qualitative features of the curves are maintained. 

In Figure 2 we present the percent error as a 
function of AVg/Vg using a resonance vc that is far 
from vt> (The 6tn mechanical resonance above was 
chosen.) The shapes of the various error curves are 
similar to those presented in Figure 1, with the 


results of the present work offering a dramatic improve- 
ment over previous approximations. 



Figure 2. Absolute value of the percent error in 
AVg/Vg versus Avg/vs, with Vc chosen to be the sixth 
mechanical resonance above 

The calculations required for the use of the 
improved formula for the dispersion are more complex 
than those required for the use of previous approxi- 
mations. An examination of the behavior of the un- 
corrected formula [Eq. (3) and Figs. 1 and 2], however, 
indicates that under a variety of conditions Eqs. (9)- 
(11) need only be applied once for any given experiment. 
One selects data corresponding to a small value of 
dispersion and uses Eqs. (9)-(ll) to compute an approxi- 
mate value of AVg/Vg. The nearly exact value for 
AVc/Vg provided by this single application of Eqs. (9)- 
(11) yields a simple multiplicative factor which can be 
applied to the mechanical resonance frequency shifts 
(Av/v) to obtain the dispersion. (This is analogous to 
the use of the (1+6) formula, with an "effective 6" 
obtained using Eqs. (9)-(ll) for a particular set of 
experimental parameters.) 

, In Figures 3 and 4 are presented the percent 
errors of the various approximations as functions of 
the transducer loading parameter 6 for two discrete 
choices of AVg/Vg.- Avg/vg large (10“^), and AVg/v^ 
small (5 x 10*®). As can be anticipated from the 
horizontal curves in Figures 1 and 2, the behavior of 
the un corrected formula and the (1+6) formula is 
independent of the size of Avg/vg. the mechanical 
resonances were chosen to be the 1st and 6th above 
in Figures 3 and 4, respectively. (The cusp-like 
behavior near 6 = 0.04 in Figure 4 is due to a change 
in sign of the error.) 
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Figure 3. Absolute value of the percent error in 
Avs/vs versus 6 = p+X-t/Ps^s Avg/vs large OO"^) 
and small (5 x 10“^ J. vc close to vt« 



Figure 4, Absolute value of the percent error In 
AVj/Vj versus 6 for Avs/vs large and small. Vr far 
from vt- 

The absolute values of the percent errors for 
the several approximations are plotted in Figures 5, 6, 
and 7 as functions of the frequency vc relative to the 
unloaded transducer resonance frequency vt (taken as 
5 MHz here). (Although smooth curves are shown, in a 
particular experiment only discrete values of Vq occur* 
corresponding to peaks of mechanical resonances. The 
resonances are spaced at 200 kHz intervals for the 
parameters used in Figures 5, 6, and 7.) Results are 


presented for small dispersion with small and large 
values of 6 1n Figure 5, and for large dispersion with 
6 ’ 0.005 In Figure 6 and with 6 = 0.2 in Figure 7. 

The results of the present work are seen to represent 
substantial improvements in accuracy for all choices 
of parameters. 



Figure 5. Absolute value of the percent error in 
AVj/Vg versus mechanical resonance frequency vr for 
AVg/vs = 5 X 10”*, 



Figure 6. Absolute value of the percent error in 
AV5/V5 versus mechanical resonance frequency vc with 
AVs/Vj large (10"^) and 6 small (0.005). 
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Figure 7. Absolute value of the percent error In 
Avc/Vj versus frequency with Avg/vs large (10“^) 
ana $ large (0.2). 


As an exaitiple of the use of the present work 
to analyze the results of an actual experiment, we 
consider magnetoelastic dispersion measurements in 
single crystal Ni.^ The specimen under study is a 
disk of length 0.1887 cm with a diameter to thickness 
ratio of about 10:1. The resulting <S for a 10.00 Wz 
quartz transducer 1s 0.0311. Crystalline axes [001], 
<111>, and [iTO], and an external magnetic field Hq 
lie in the plane of the disk. We define 6 q as the 
angle between Hq and the [001] direction. Transverse 
ultrasonic waves of frequency '^8.9 IWz were propagated 
along the [110] axis perpendicular to the plane of the 
disk. The frequency vc of a particular mechanical 
resonance was measured as a function of Ho and 0 q. 

The ultrasonic dispersion Avg/vg was calculated using 
6 q • 90°, Hq = 10 kOe as the uncoupled magnetic field 
orientation. In Figure 8, the data have been reduced 
using the un corrected formula and the present work. 
Even on this linear plot, the improvement provided by 
the present work is evident. 
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Figure 8. Ultrasonic dispersion (Avg/vg) in bulk 
single crystal Ni versus magnetic field Hq for 
6g = 0° and 00 = 90°. Dispersion determined using the 
uncorrected formula and the present work. 
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